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ORIGINAL PAGE IS 

OF POOR QUALITY SUMMARY 

Closed form solutions exist In the literature for propagating plane har- 
monic waves through infinitely long uniformly thick plates which are immersed 
in fluids. These solutions fail to account for the edge effects present in 
finite length plates and address only plates of uniform thickness. This paper 
presents an analysis using finite element techniques which addresses the prop- 
agation of a uniform incident pressure wave through a finite diameter axisymme- 
tric tapered plate immersed in a fluid. The approach utilized in developing a 
finite element solution to this problem is based upon a technique for axi symme- 
tric fluid structure interaction problems presented in the Tenth NASTRAN Con- 
ference (ref. 1). 

The problem addressed in this paper is that of a 10 inch diameter axi sym- 
metric fixed plate totally immersed in a fluid. The plate increases in thick- 
ness from approximately 0.01 inches thick at the center to 0.421 inches thick 
at a radius of 5 inches. Against each face of the tapered plate a cylindrical 
fluid volume has been represented extending five wavelengths off the plate in 
the axial directon. The outer boundary of the fluid and plate regions have 
been represented as a rigid encasement cylinder as was nearly the case in the 
physical problem. 

The primary objective of the analysis is to determine the form of the 
transmitted pressure distribution on the downstream side of the plate. A 
secondary purpose is to demonstrate the feasibility of employing the ref. 1 
fluid-structure interaction technique on a large D.O.F. problem, and in fact 
running this large problem uncovered a basic error in the NASTRAN coding for 32 
bit machines ( e . g . , VAX or IBM 360, 370, or 3033) when employing MPC or DMIG 
cards. These errors could impact any rigid format. 

The finite element model used consists of 7455 nodes comprising 7191 axi- 
symmetric ring elements. Both the fluid and the plate structure have been 
represented using axisymmetric CTRAPAX and CTRIAAX elements. The fluid (sea- 
water) elements and structure (plate) elements are linked together along their 
boundary by utilizing a duplicate set of nodes and a direct matrix insertion 
capability within NASTRAN. This enables the incorporation of the stiffness and 
mass coupling terms in the stiffness and mass matrices respectively. 

The plate is loaded with a uniform harmonic pressure wave (f = 100,000 Hz) 
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which propagates axially through the Incident side fluid, Into the plate, and 
then out Into the transmitted side fluid. 

Presented results consist of plots of pressure along both the Incident and 
transmitted surfaces of the tapered plate. In addition, pressure and piiase 
contours throughout the modeled fluids are presented. 


SYMBOLS 

Values are given in both SI and U.S. Customary Units. The measurements 
and calculations were made in U.S. Customary Units. 

AA ■ associated area (eq. (1)) 

c ■ value of term entered on diagonal of damping matrix 

c f * sound wave velocity in fluid 

f * frequency of propagating wave 

kf = bulk modulus of the acoustic fluid 

ry * normal to a surface In radial direction 

ry » normal to a surface In axial direction 

U 2 * displacement In axial direction 

Pf * mass density of acoustic fluid 


INTRODUCTION 

Sternberg, et al., (ref. 2) have experimentally demonstrated the frequency 
Independent characteristics of acoustic filter plates with regards to beamform- 
ing characteristics. Mathematical approaches have been undertaken to analyti- 
cally account for this phenomena with vary I t degrees of success. The objec- 
tive of this project was to utilize the finite element method to model the 
transmission of a plane harmonic wave through a filter plate of variable thick- 
ness Immersed in sea water. This filter plate effect, If used In conjunction 
with an acoustic lens, develops beams which are Independent of frequency and of 
constant width. A typical configuration Is shown in Figure 1 and has been 
experimentally tested. 

The analysis presented accounts for the rigid containment cylinder, the 
tapered circular steel filter plate, the transmission side fluid Inside the 
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containment cylinder and the Incident side fluid which theoretically extends to 
Infinity. The rubber acoustic lens discussed briefly above has been excluded 
from this analysis in order to Isolate the effects of the filter plate. 

Since the actual structure Is axlsymmetrl.* the finite element model uti- 
lized Is also axl symmetric. The filter plate has been analyzed by using the 
displacement method within NASTRAN. The fluid has been modeled In accordance 
with an axl symmetric pressure analogy (ref. The driving force consists of 
a plane harmonic wave at 100 kHz at normal Incidence to the steel plate. Pres- 
sure and phase distributions throughout the Incident and transmitted fluid 
volumes are to be computed. To enhance the presentation of the rather exten- 
sive pressure and phase data, pressure and phase contours throughout the 
Internal and external fluid have been presented. In addition, plots of pres- 
sure along the Incident and transmitted surfaces of the filter plate are 
presented. 


Model Development 

As detailed In ref. 1 the development of a finite element model for any 
axl symmetric fluid structure Interaction problem requires several steps. For 
large degree of freedom problems model sizing, element sizing, and element 
selection are reasonably straightforward, however, they must be addressed cau- 
tiously In order to avoid both technical mlsjudgments and excessive computer 
cost. 


I. Sizing . The particular piece of hardware modeled Is well defined. As 
shown In figure 2 the encasement cylinder Is 5" (12.7 cm) In length and 0.5" 
(1.27 cm) thick. The filter plate Is 10" (25.4 cm) In diameter and tapers from 
.08" (0.203 cm) at the center to 0.421" (1.06 cm) at a 4" (10.16 cm) radius. 
The location of the driving surface Is 2.5" (6.35 cm) away from the filter 
plate surface and at the end of the containment cylinder. The area between the 
filter plate and driving surface Is fluid filled. On the downstream side the 
modeled fluid extends to the end of the containment cylinder. The fluid within 
the containment cylinder extends 4 wavelengths off the plate. 

The Important aspect relating to element sizing for steady state fluid 
structure Interaction problems Is that the distribution of nodes throughout the 
model be such that plane harmonic waves at .lie frequency of Interest can ade- 
quately be described. As developed In ref. 3 adequate representation of the 
pressure waves Is ensured by using ten nodes per wavelength or In this case 
where f 3 100 kHz the node spacing becomes 0.06" (0.152 cm). This Is based 
upon a wave propagation velocity In seawater of 60000 In/sec (152400 cm/sec). 
With node spaclngs this small and a model geometry as described above, the num- 
ber of degrees of freedom approaches the limits of COSMIC NASTRAN on the VAX 
11/780. 

Since an axlsymmetrlc pressure analogy solution has been Invoked for the 
fluid elements, the selection of mating element type Is limited. Several axl- 
symmetrlc elements exist within COSMIC NASTRAN, however, previous success with 
CTRAPAX and CTRIAAX elements determined selection. These elements have been 
used to represent both the fluid and filter plate. As will be discussed In the 
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next section, the containment cylinder and driving surfaces are represented as 
boundary or load conditions. The finite element model developed according to 
the above discussion is shown in Figure 3. 


II. Boundary Condi ti o n . The boundary condition addressed next Is the 
coupling of the fluid elements (modeled within an axi symmetric pressure 
analogy) to the structure elements (modeled within the displacement method). 
This coupling when using axisymmetric elements requires duplicate nodes along 
any fluid structure boundary. These duplicate nodes are then utilized In 
formulating the off diagonal coupling terms within both the mass and stiffness 
matrices for each of the elements involved. The details of the coupling pr- 
ocedure are given in ref. 1 and are not repeated here. 

w w 

The mass terms are given as, - A n r and - A k (1) 

and stiffness terms are given as, AA and AA n 

The physical nature of the filter plate hardware enables representation of 
the cylindrical containment cylinder as a rigid outer fluid boundary condition, 
and is expressed as 



where pis the fluid pressure. 

This boundary condition is implemented within the pressure analogy by 
freeing all constraints along the containment cylinder/fluid Interface (i.e., 
no surface tractions). 

The external fluid boundary is terminated with fluid dampers in sucn a way 
that the fluid appears to extend to mathematical infinity. Since NASTRAN does 
not contain any axisymmetric dashpot elements, this effect was developed using 
a direct matrix insertion technique similar to that required along the fluid/ 
structure boundary. To each node along the fluid boundary the magnitude of the 
lumped damping term to be inserted in the damping matrix is given as 

AA b k f 


As with the stiffness and mass terms described previously, the location of the 
damping term within the damping matrix is developed in ref. 1. 

III. Loading . Plane harmonic pressure waves will be the predominant load 
type seen by the filter plate. The plane harmonic wave utilized In this prob- 
lem is at a frequency of 100 kHz and has a magnitude of 1 psl (6895 P ). Ti.ls 
wave emanates from the outer fluid boundary parallel to cfce flat metal surface 
and Impinges on the filter plate at normal incidence. The boundary conditions 
and loading scheme are shown in Figure 3. 
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DISCUSSION 


Initial attempts to execute NASTRAN and develop solutions for the pressure 
field resulting from the tapered filter plate Mere unsuccessful. Several out 
of the ordinary problems Mere encountered Mhile making these Initial checkout 
computer runs. The first of these relates to the use of the DMIG card or W>C 
card option Mlth larger problems on a 32 bit machine (In this case the VAX 
11/780). 

After considerable time and expense. It Mas discovered that the problem 
arises Mhen a model has a large numbegpof nodes such that the combined roM and 
column Internal node number exceeds (the 32 bit Mord can no longer repre- 
sent an unsigned Integer of this size). As a result Mhen NASTRAN attempts to 
order the list of roM *nd column ID's In a monitonlcally Increasing order, the 
large Internal node numbers actually appear as smaller negative numbers because 
the sorting routine processes the 32nd bit as a plus or minus sign. As a 
result the DMIG card coupling terms Mere arbitrarily associated Mlth a node 
elseMhere In the matrix thus voiding the coupling effect. Coding changes Mere 
suggested by RPK Corporation and ultimately corrected this error. The coding 
change alters the Internal node sorting routine In such a May that the 32nd bit 
Is read as part of the node number and not as a plus {+) or minus (-) sign. 
Subsequent runs deinonstrated proper placement of the DMIG card coupling terms. 

The second, less serious p" ! »n encountered relates to the mass distribu- 
tion for elements Mhich are very ciose to the centerline. 

NASTRAN calculates the lumped massof any axisymmetric element by first determin- 
ing the centroid of the elements nodes. This radius is then used to determine 
the circumference of the ring element. KnoMlng the length and height of the 
element then enables a determination of the area and therefore mass of the 
entire ring element. This mass is allocated equally to each of the four nodes 
comprising the elements connectivity. For large radii this allocation sc erne 
is very accurate. iiOMOver, Mhen the radius to che centroid of the element is 
less than five times the node spacing significant errors result (as much as 
40%). This error has been corrected by using the elements outer and Inner 
radii to calculate ring areas relative to the centroidal radius. This divides 
each ring element into an outer annul us and an inner annulus. The mass of the 
outer annulus is allocated to the tMO outer annulus nodes. The mass of the 
inner annulus is allocated to the tMO inner nodes. The FORTRAN coding changes 
made to the TRAPAD and TRIAAD subroutine in NASTRAN are shOMn in Figures 6 and 
5. The lines that have changes are mainly in the lines Mi thou t the column 
73-80 card number designators. 


RESULTS 

The NASTRAN solution to a steady state fluid structure Interaction problem 
consists of the pressure and phase distribution throughout the fluid. For very 
large problems this data is ^'mbersome to Interpret, consequently tMO data pre- 
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sentatlons have been developed. The first consists of a plot of the pressure 
along the upstream and downstream sides of the tapered metal plate. This data 
shown in Figures 6 anc 7 show several areas of higher pressure transmission 
through the tapered plate. The major peak to peak fluxuations along the plate 
surface in the pressure plots correlate with the submerged plate flexure wave- 
lengths existing In a constant- thickness submerged plate whose thickness is the 
average type thickness. 

A finer mesh distribution radially along the centerline is needed in order 
to properly model the stiffness in the plate and thus pick up the transmitted 
pressure. This is why the transmitted and incident pressures are too large at 
the centerline. Subsequent runs with a fluid plate (field solution is known) 
demonstrated that this error near the centerline was localized to two or three 
elements. 

A second and, in this case, more useful presentation of this data <s a 
color graphics format of pressure and phase contours. Pressure contours are 
shown in Figure 8. A noticeable effect depicted by these pressure contours is 
the redirection or focusing of the pressure towards the centerline. Equally 
important, the appearance of annular regions of increased energy transmission 
are evident. There are regions where more energy transmits through the plate 
than in similar neighboring regions. The phase contour shown in Figure 9 is 
evidence of the transmission of an actual plane wave from the driving boundary 
through the incident side fluid. The bands of constant and repeating phase 
angles verify this plane wave. The discontinuities in the phase contours in 
the downstream fluid correlate with certain redirected pressure contours evi- 
dent in Figure 8. These acoustic characteristics of filter plates, while 
experimentally verified, have not been readily predictable. The technique 
utilized here establishes a reasonable prediction capability. The development 
of a predictive capability for both the redirection of the pressure field and 
the regions of higher energy transmission were objectives of this effort. 
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COUPLE i FLUID AXIS 0 F 



FIGURE 3. A FINITE ELEMENT MODEL OF THE TAPERED FILTER PLATE, 
SURROUNDING FLUID, AND STRUCTURAL BOUNDARIES 
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FIGURE E FORTRAN CODING CHANGES IN NASTRAN SUBROUTINE TRIAAD FOR CTRIAAX ELEMENTS 
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OF THE TAPERED FILTER PLATE 


00000 * 



FIGURE 7 THE PRESSURE DISTRIBUTION AS A FUNCTION OF RADIUS ALONG THE TRANSMITTED 
FACE OF THE TAPERED FILTER PLATE 







